Spin Excitation Assisted by Non-Softening Phonon 
for Spin-Peierls Model 
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We study spin dynamics of a spin-Peierls chain with nearest-neighbor and next-nearest-neighbor 
Heisenberg spin exchange interactions together with a gapped and dispersionless phonon. The dy- 
namical spin correlation function and phonon excitation spectrum are calculated at zero temperature 
by using dynamical density-matrix renormalization-group method. We find a new spin excitation 
assisted by non-softening phonon. The excitation is located above phonon in energy and shows a dis- 
persive feature with strong intensity near the momentum it. The phonon excitation spectrum is also 
influenced by the spin-phonon interaction. We discuss the possibility of observing the spin-phonon 
coupled features in inorganic spin-Peierls compound CuGeC>3. 



I. INTRODUCTION 

One-dimensional (ID) quantum spin system coupled 
with lattice degree of freedom has been extensively stud- 
ied experimentally and theoretically, since the systems 
provide a playground of spin-Peierls transition. Conven- 
tional spin-Peierls compounds like organic materials ex- 
hibit the transition with spin dimcrization and lattice al- 
ternation accompanied by soft- phonon mode. Theo- 
retically the spin-Peierls transition can be derived under 
the presence of soft-mode phonon by using the random 
phase approximation approach in the adiabatic limit. Q 
The discovery of CuGe03, [|| however, has casted a prob- 
lem on the conventional mechanism of the spin-Peierls 
transition, since the soft-phonon mode associated to lat- 
tice alternation has never been found in CuGe03 so 
far. [3-0] 

A recent study has shown a theoretical explanation of 
spin-Peierls instability in the antiadiabatic limit. Q The 
spin-Peierls Hamiltonian can be mapped approximately 
to spin-i J1-J2 Heiscnberg Hamiltonian by the flow equa- 
tion in the antiadiabatic limit, where J\ and J2 is nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) ex- 
change interactions, respectively. The spin-Peierls tran- 
sition, therefore, corresponds to the Kosterlitz-Thouless 
(KT) transition in the J1-J2 Heisenberg model obtained 
by changing the ratio of J\ and J2. The KT transi- 
tion is order-disorder type quantum phase transition, and 
has a critical exchange ratio, a c = 0.241167, between 
spin-liquid and dimer phases at the ground state. [9l-fli1| 
Namely, for the exchange ratio J2/J1 = ct < a c , the 
ground state belongs to the spin-liquid phase, but for 
a > a c , it is in the dimer phase. The renormalization of 
phonon degree of freedom by the flow equation changes 
the ratio a: a increases as spin-phonon coupling g in- 
creases. Thus, there occurs the KT transition without 
soft phonons in the antiadiabatic limit for large spin- 
phonon coupling g > g c = 0.8. [l2T - fl8T | This means that 
no soft phonons exist below the critical temperature. [|[ 
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The spin-Peierls transition was experimentally ob- 
served for the organic compounds, tctrathiafulvalene 
(TTF), tetracyanoquinodimethane (TCNQ) series in 
the 1970s, P, 13, flW22l| and the inorganic compounds 
CuGcC>3 in 1993. [J] Experimental studies have shown 
that phonon frequency associated with lattice distor- 
tion in CuGeOs is higher than that in the organic com- 
pounds. In (TTF)CuS4C 4 (CF 3 )4, a soft-phonon fre- 
quency is estimated to be w p h = 1.4 me V = 16 K, and 
spin-Peierls gap is A = 21 K. [l], [l9| On the other hand, 
for CuGe0 3 , a dispersive phonon frequency related to 
lattice distortion is cjph = 6.8 THz = 330 K, a disper- 
sionless phonon is w p h = 3.2 THz = 150 K, and the gap 
is A = 2 meV = 23 K. @, A condition uj ph < A 

for the conventional spin-Peierls mechanism is fulfilled in 
the organic compounds, but not in CuGeOs. Therefore, 
the spin-Peierls transition in CuGeOs can be explained 
by the theory in the antiadiabatic limit. 

CuGeOs has an antiferromagnetic NN interaction J\ 
and an antiferromagnetic NNN interaction J2 that in- 
duces frustration. Experimental and theoretical studies 
have shown that J x w 100 K and a = 0.36. 0-0] The 
instability of the spin-Peierls phase due to J2 has been 
discussed by WeiBe [26| . 

Spin excitations in CuGeOs have been observed by in- 
elastic neutron scattering. |27H2g | The separation of spin 
and phonon excitations in the experiments has not been 
complete, although some of phonon excitations have been 
identified by the same technique. p3l - [25| A recent devel- 
opment of polarized neutron scattering may resolve this 
problem in the near future. To give theoretical supports 
for inelastic polarized neutron scattering experiment, we 
investigate spin excitation of the spin-Peierls model with 
the non-softening phonon. Our chief aim in the present 
paper is to clarify the effect of non-softening phonon on 
spin excitations in CuGeOs. 

We use the dynamical density-matrix renormalization- 
group (D-DMRG) method to calculate dynamical 
spin correlation function and phonon excitation spec- 
trum. This method is a dynamical version of DMRG 
method [3(| presented by Jeckelmann [3lT |. 

In this paper, we perform D-DMRG calculation of the 



2 



dynamical spin correlation function for both the spin- 
Peierls model and its effective spin model (J1-J2 model). 
Treating phonons in the spin-Pcicrls model as quantum 
objects, we find a new spin excitation assisted by non- 
softening phonon. There is no corresponding structure 
in the effective J\-Ji model. The new structure shows 
a dispersive feature with strong intensity near the mo- 
mentum 7r, and it is located above the phonon energy. 
The new structure is explained by using particle-hole ex- 
citation assisted by phonon in the Su-Schricffcr-Hccgcr 
model onto which the spin-Peierls model with XY-type 
spin is mapped. The phonon excitation spectrum is also 
influenced by the spin-phonon interaction. Main peaks 
of phonon near the momentum tt show a new tail struc- 
ture whose energy range is the same as the new spin 
excitation. In addition, spin-assisted phonon structures 
appear at low-energy region. These features induced by 
spin-phonon interaction are expected to be observed in 
inelastic neutron scattering experiments in the near fu- 
ture. 

This paper is organized as follows. In Sec. ITTJ we show 
a model Hamiltonian of the frustrated spin-Peierls chain 
with Einstein phonon and introduce a rcnormalized J±- 
J2 Hamiltonian. We explain the numerical method, D- 
DMRG, in Sec. IIIII In preparation for the spin-Peierls 
model, we demonstrate dynamical spin correlation func- 
tion in the J1-J2 model in Sec. IIV Al Dynamical spin 
correlation function and phonon excitation spectrum of 
the spin-Peierls model are presented in ElIVBl Summary 
of the present paper is given in Sec. [V] 



II. MODEL 

We consider the following model Hamiltonian that de- 
scribes ID frustrated spin-Peierls chain. 



with 



tL — H s + Hp + %sp 
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H tv =g^2(x j -x j+1 )hf\A), 



and 



hP(A) = S j -S. 
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(1) 

(2) 

(3) 
(4) 

(5) 



where Sj is the 7-site spin operator for S = |, x 3 is 
the j-site displacement of the coordinate with respect to 
equilibrium position, and pj is the conjugate momentum 
operator. J is bare NN exchange interaction, and ao is 
the ratio of bare NN and NNN interactions. In this work, 
we consider only antifcrromagnctic coupling for the NN 



and NNN interactions. M and K are the effective mass 
and the elastic coupling constant, respectively, g is the 
spin-phonon coupling constant. A is the Ising anisotropic 
parameter. 

After the second quantization of lattice degree of free- 
dom, the phonon part of the Hamiltonian ([3]) and the 
spin-phonon coupling term are rewritten by 



and 



AJ 
~2 



J 3 + l 



(0) 



h ( - ] (A), (7) 



respectively, where bj and bj are the creation and annihi- 
lation operator of j-site phonon, respectively. We em ploy 
the Einstein phonon with frequency luq = y K/M. The 
dimensionlcss spin-phonon coupling constant A is defined 
as A = {g/ J)y/2/MujQ. We use the natural unit, h = 1. 

The spin-Pcicrls model without the NNN interaction in 
the XY limit (ao = and A = 1), that we call XY spin- 
Peierls model, can be mapped to the Su-Schrieffer-Heeger 
(SSH) model by the Jordan- Wigner transformation, 



H e = -J^cos(fc)4cfc 



and 



H cp = iXJ 

k,l 



sin(fc) - sin(Z)] (b k -i 



(8) 



(9) 



where b q and are the momentum representation of the 

phonon creation and annihilation operators, ct and Ck 
are the momentum representation of the spinless-chargc 
creation and annihilation operators. 

An effective Hamiltonian with Heisenberg spin ex- 
change interactions (A = 0) after renormalization of 
phonon degree of freedom by the flow equation is given 
by 1 



N 



N 



H 



J1J2 



(10) 



with 



Ji = J 
J 2 = J 



A 2 J 3(l-a )A 2 J 2 



a 



A 2 J (3-5a )A 2 J 2 



8w 2 



+ 0(A 3 ), (11) 
+ 0(A 3 ), (12) 



where \J /loq and J/wo are small parameters for expan- 
sion and 0(A 3 ) is the third or higher order terms of A. 
Equation (fT0|) is a 3\-Ji model with effective NN and 
NNN interactions. In this model, with increasing the ra- 
tio a = J2/J1 from zero, the KT transition from spin 
liquid phase to dimer phase occurs at a c . [9MTT| In the 
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spin-liquid phase where a < a c , a gapless excitation ap- 
pears. On the other hand, in the dimer phase where 
a > a c , there is a spin gap between the ground state 
and the lowest triplet excitation, accompanied by spon- 
taneously dimerization of spin pairs. The quantum phase 
transition does not require any soft phonon. 



III. METHOD 

To examine dynamical behavior for the spin-Peierls 
model, we calculate two quantities, dynamical spin cor- 
relation function and phonon excitation spectrum. 

The dynamical spin correlation function at zero tem- 
perature is given by 



Xs (fc, w ) = S-L(0|S z (fc)- 
ttL lo 



1 



+ e - % + 17 



S z (k)\0), (13) 



where L is the system size, 7 is the damping factor with 
small positive number, and |0) is the ground state. S z (k) 
is the z component of the spin operator at momentum k. 
We define phonon excitation spectrum as 

Xp (q,u)) = -%±(0\b(q)— 1 -——b\q)\0), (14) 

ttL oj + €q — H + vy 

where b(q) and b^(q) are the phonon annihilation and cre- 
ation operators at momentum q, respectively In eq. (|14p . 
we consider a pair of b(q) and b'(q). In the phonon 
Green's function, there are additionally three pairs of 
b^(q) and b(q), b(q) and b(q), and b^(q) and b^(q). If the 
spin-phonon coupling is small, there are few phonons in 
the ground state. In such a case, the contribution from 
the three pairs is expected to be small. Therefore, we 
consider only the case of eq. (|T4|) to represent the phonon 
excitation spectrum. 

For strongly-correlated ID systems, it is well-known 
that the DMRG method can pro vide a good numerical 
solution of the ground state. [3(| The dynamically ex- 
tended version of DMRG, D-DMRG, is also suitable for 
obtaining dynamical properties at zero temperature. [3l[ 
Since the spin-phonon interaction breaks the mir- 
ror symmetry, we consider two reduced density matri- 
ces for both the system and the environment blocks in 
the DMRG process. We use the infinite-size algorithm 
of DMRG [30[ for systems with open boundary condition 
(OBC). In OBC, the Fourier transform of the spin and 
phonon operators read 



s z (k) = J J2 S 1MM 



(15) 



L 



b (q) = y^rl]^ sin (i?). 



(16) 



with momenta k and q given by mr/(L + 1), (n 
1,2,-.. ,L). 



For the D-DMRG method, we use three target 
states: |0), S z (k)\0), and [u + e - H + i-y]' 1 S z (k)\0) 
for eq. (TBI, and |0), 6 + (<?)|0>, and 
[u + eo-H + iil^tfiq^O) for eq. KHJ). In the D- 
DMRG procedure, we use a modified version of the 
conjugate gradient method to calculate the correction 
vector, [ui + e - U + ry] -1 6\Q). 



IV. RESULTS 

In this section, we first show calculated results of dy- 
namical spin correlation function for the effective J1-J2 
model (|10[) deduced from the spin-Peierls Hamiltonian 
(P). Secondly, treating phonons quantum-mechanically, 
we calculate the dynamical spin correlation function of 
the spin-Peierls model ([T]). New structures originated 
from the quantum phonons arc identified by making a 
comparison with the J\-Ji model. Finally the effect of 
spin-phonon coupling on phonon excitation spectrum is 
examined. 

In order to simulate CuGe03, we take the phonon en- 
ergy ujo = 3 J, [23M25I ] except in the case explicitly 
provided. The value of spin-phonon coupling is not clear 
for CuGe03. Therefore, we take A satisfying AJ/wq < 1. 



A. Effective J1-J2 model 

The dynamical spin correlation function for a 64-site 
J1-J2 Hciscnbcrg chain is shown in Fig. [TJ The damp- 
ing factor 7 in eq. (fT3|) is taken to be 0.1 J. We note 
that a preliminary result has been reported in rcf. 32). 
In Fig. [lfa), we take cxq = and A = 1.12, resulting in 
J2/J1 =0.1 from eqs. (HI]) and (jT2j). The ratio is below 
a c . This means that Fig. [2(a) represents the dynamical 
spin correlation in a spin-liquid phase. The distribution 
of spectral weight is similar to the exact results of the 
ID Hciscnbcrg model (J2/J1 = 0) where spectral weight 
consists of the des Cloizeaux-Pearson mode at the lowest- 
energy branch and multi-spinon continuum. [33T - I351 ] The 
inset in Fig.QJa) shows the system-size dependence of the 
position of the des Cloizeaux-Pearson mode at the small- 
est momentum fe = 7r/(£ + l). A fitting function gives 
nearly zero excitation energy at l/L — > as expected in 
the spin-liquid phase. 

The spin correlation function in a dimer phase is shown 
in Fig.QTb), where a a = 0.36, A = 1.06, and thus 

J2/J1 = 0.4. We can find two characteristics as com- 
pared with Fig.rjja), i.e., strong intensity around k = tt/2 
(ui/Jr^j 1) and a peak structure at the upper edge of the 
spinon continuum around k = ir (cj/J ~ 2). These fea- 
tures in the dimer phase have been reported in a previ- 
ous study for small systems up to 16 sites under periodic 
boundary condition. 0, H?} The inset shows the pres- 
ence of spin gap in the thermodynamic limit. The gap 
magnitude (~ 0.08 J) is similar to a previous report eval- 
uated from DMRG calculations. 381 These results also 
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4 single-site dimension of spin-phonon coupled system is 

several times larger than that of pure spin system, result- 
ing in the requirement of huge computational resources. 
SCF ^ e cnec ' teo - t ne convergence of the calculation in terms 
of two parameters, i.e., the DMRG truncation number, 

2 to, and the maximal phonon number per site, n p . We 

found that good convergence not only for the ground 
state but also for spectral weight at high-energy region 
can be achieved for to = 100 and n p = 2. In this paper, 
we use the broadening factor for the dynamical quantities 

7 = 0.1.7. 
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FIG. 1. (Color online) Intensity map of dynamical spin cor- 
relation function (SCF) in a 64-site J1-J2 chain under two 
conditions: (a) J2/J1 = 0.1 {loo/ J — 3, Qo = 0, A = 1.12), 
(b) J2/J1 = 0.4 (ojo/J = 3,a = 0.36, A = 1.06). The inset 
in both panels shows system-size dependence of the energy of 
peak position at the smallest momentum k = n/(L+ 1). Blue 
squares represent calculated peak positions, and red lines de- 
note fitting by a function with a power: /(1/L) = a(l/L) b + c. 
a = 5.63, b = 0.90, and c = -0.01 for (a), and a = 7.43, 
b = 1.14, and c = 0.07 for (b). The gap in (b) remains finite 
in the thermodynamical limit 1/L — > 0. The truncation num- 
ber of DMRG and the broadening factor of Lorentzian are set 
to be m — 100 and 7 = 0.1J, respectively. 

confirm the validity of our D-DMRG calculation. 

There is a breakpoint near k = -k/2 in the lowest- 
energy branch in Fig. |TJb). This is an artifact by finite- 
size effect under OBC, since momenta near k = ir/2 have 
a substantial contribution from the edges of the system 
as expected from eq. (|15p . This, however, does not occur 
for momenta close to k = and close to k = n, for 
which sin(jfc) in eq. (| 15|) goes to zero with approaching 
the edges of the system. 



B. Spin-Peierls model 

We calculate dynamical spin correlation function and 
phonon excitation spectrum for a 16-site frustrated spin- 
Peierls chain in this section. The system size is smaller 
than that for the J1-J2 model (64 sites). This is because 



1. Dynamical spin correlation function 

Figure [Ha) shows the case without spin-phonon cou- 
pling, i.e., only H s term in eq. ^ with ao = 0.36. This 
is nothing but a 16-site result of the J1-J2 model, and 
the spectral behavior is consistent with the 64-site result 
with spin gap mentioned in QYV Al 

The introduction of A changes the spectrum as shown 
in Fig. [U^b), where XJ/uj = 0.5. The most striking 
change is the appearance of high-energy spectral weight 
around u> ~ 4J (see red solid lines and intensity map 
shown in the inset). The intensity is strongest at k = 7r 
and widens toward low-energy side with decreasing k fol- 
lowed by an energy minimum at k = n/2. The energy 
position is higher than the dispersionless phonon located 
at wo = 3J. From this result, it is expected that addi- 
tional spin excitations originating from the spin-phonon 
coupling exist above the energy of the phonons. Note 
that the intensity of the new structure increases with in- 
creasing A. 

Comparing red solid lines in Fig. EJb) with those in 
Fig.[^a), we find a change of spectral distribution along 
the lowest energy branch: the spectral intensity at k = ti 
is suppressed and the weight transfers toward k = n/2. 
This change is similar to that caused by increasing J2/ J\. 
This is reasonable, since the effective value of J2/J1 eval- 
uated from the second-order expression of A in cqs. (|lip 
and (fT2|) is a = J2/J1 = 0.44, which is larger than the 
bare value cto = 0.36. For comparison, the spin cor- 
relation function for the J1-J2 model with a = 0.44 is 
plotted in Fig. [2fb) as blue dotted lines. The lowest- 
energy branch follows that of the spin-Peierls model. 
This tempts us to justify the use of effective J1-J2 model. 
However, we can find a qualitative difference between 
the J\-J-i model and the spin-Peierls model: the upper 
edge of multi-spinon excitation (~ 2.5 J near k = n) in- 
creases in the spin-Peierls model while decreases in the 
J1-J2 model as compared with that of Fig.[2fa). A simple 
explanation of the difference would be that the second- 
order contributions in eqs. (fTT|) and ([T2|) are not enough 
for the complete description of the spin-Peierls model 
and higher-order terms contributes significantly for the 
present parameter set. 

Figure [^c) shows the case of ao = but with the same 
A as Fig.[^b). As is the case of Fig.[^b), there appears 
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FIG. 2. (Color online) Dynamical spin correlation function 
(SCF) in a 16-site spin-Peierls chain. The phonon energy 
ojo = 3J. (a) ao = 0.36, AJ/wo = 0, (b) ao = 0.36, XJ/loq = 
0.5, and (c) ao = 0, XJ/uo = 0.5. The blue dotted lines 
in (b) and (c) show the dynamical spin correlation function 
of the effective J\-Ji model, where Ji and J2 are evaluated 
from eqs. (fTT) and (fT2|) . respectively: J2/J1 = 0.44 in (b) 
and J2/J1 = 0.17 in (c). The inset in each panel shows the 
intensity map of the correlation function. The truncation 
number of DMRG is set on m — 100 and maximal phonon 
number is two. The broadening factor is set to be 7 = 0.1J. 




FIG. 3. (Color online) Dynamical spin correlation function 
in a 16-site spin-Peierls chain with ao = 0.36 and uq/J = 
1.5. The red lines represent the case of AJ/wo = 0.5, while 
the blue dotted lines represent the case without the coupling. 
The inset shows intensity map of the correlation function for 
AJ/cjq = 0.5. 



a high-energy dispersive structure with small intensity 
around u> = 4J induced by the spin-phonon coupling in 
the spin-Peierls model. The upper edge of multi-spinon 
excitations (u> < 3J) in the spin-Peierls model is larger 
than that of the effective J1-J2 model. This is again the 
same as Fig. EJb), indicating insufficient mapping of the 
spin-Peierls model onto the J1-J2 model for the present 
parameter set. We find small spectral weights below the 
des Cloizeaux-Pearson mode. The weights are caused by 
finite-size effect and decrease with increasing system size. 

In order to examine the case where phonons are in- 
side the multi-spinon continuum, we take ujq/J = 1.5 
(refs. 23-25) and show the spin correlation function for 
XJ/luq = 0.5 (red solid lines) in Fig. Making a com- 
parison with the case without the coupling (blue dot- 
ted lines), we find that phonon- induced spin excitations 
appear at u ~ 2. 5 J with a dispersive structure show- 
ing a minimum at k = tt/2. We note that, although 
the coupling constant A = 0.75 is smaller than the cases 
of Fig. [^b) (A = 1.5), the spectral intensity relative to 
multi-spinon continuum is comparable to Fig.[2^b). This 
probably comes from the enhancement of hybridization 
between spin and phonon due to overlapping of their en- 
ergy scale. 

Since the phonon-induccd spin excitation appears in 
broad energy range, it is expected to be coupled with the 
multi-spinon continuum. However, in the dimer phase 
of the spin-Peierls models, there is not only the contin- 
uum but also the state with large weight at the lower 
edge of the continuum. In order to make clear how the 
continuum contributes to the phonon-induced spin exci- 
tation, we examine the XY spin-Peierls model (ao = 
and A = 1 in eqs. @ and ((?])), where the XY spin chain 
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FIG. 4. (Color online) Dynamical spin correlation function in 
a 16-site XY spin-Peierls chain (or dynamical charge-charge 
function in a 16-site Su-Schrieffer-Heeger chain). The inset 
shows intensity map of the correlation function for AJ/wo = 
0.5 and luq/J = 3.0. 



has only continuum excitations of spinons. [37| Figure 0] 
shows the dynamical spin correlation function for the XY 
spin-Peierls model with A J/wo = 0.5 and ojq/J — 3.0. We 
can see the same structure induced by the spin-phonon 
coupling in the high energy region as the case of the spin- 
Pcierls model. Therefore, we conclude that the phonon- 
induced spin excitation couples with the continuum of 
spinons. 

Moreover, we investigate the A and loq dependence 
of the phonon-induced spin excitation in the XY spin- 
Pcierls model. We show A dependence of the dynami- 
cal spin correlation function near k = tt with XJ/ujq = 
0, 1/6, 1/3, and 1/2 for fixed w /J = 3.0 in Fig. E (a), 
and cjo dependence with uq/J = 2.0,3.0, and 4.0 for 
fixed AJ/cjo = 0.5 in Fig. [5] (b). We find that the en- 
ergy position of the excitation does not depend on the 
spin-phonon coupling A, but the integral of the excita- 
tion depends on A with a power-law behavior as seen in 
the inset of Fig. [5fa). In Fig. [IJb), we can see that the 
excitation is situated between ojq and luq + 2 J at k = tt. 

The behaviors in the dynamical spin correlation func- 
tion in the XY spin-Peierls model can be explained by 
the SSH model in eqs. <JSj> and J9]). We assume that 
phonon creation and annihilation operators never have 
any finite expectation values in the ground state. This 
assumption is supported by the well-known fact that the 
ground state is not in the dimer phase with the bond 
alternation for small coupling constant A. [4l[ The spin 
correlation function is rewritten by the spinless charge- 
charge correlation function. We obtain the imaginary 
part of the spinless charge-charge correlation function, 
X c (fc,w), within the second-order perturbation in terms 
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FIG. 5. (Color online) Dynamical spin correlation function 
in a 16-site XY spin-Peierls chain at k = 16n/17 = n. (a) 
the A dependence. The inset shows the A dependence of the 
integrated weight in the energy range of 3 J and 6 J for uio/ J = 
3.0. The blue crosses represent the integrated weight obtained 
by subtracting the integrated weight at A = 0, and the red line 
denotes the power-law fitting, a(\J/u)o) b with a = 4.4 x 10~ 2 
and b = 1.8. (b) Wo dependence for \J/luq = 0.5. The inset 
represents x£ (7t,cj) + Xc 2i> ' (w, <^) obtained by eqs. (1191) and 
(Ell). 



of A: 

Xc (fc, w) = xi 0) (k, w) + x? a) (k, u) + xi 2h) (M) + 0(A 4 ) 

(17) 

with 

Xi 0) (fc,^) = J2 6 ^ - tl + k + ^n-ei)e(e l+k ) (18) 
l 

x^froo = f^y v (^ +M+fc+g)x 

Lq 



x? b) (fc,^) 



x 6(-ei)6(t l+k )0(e l+k+q ) 
x 8{u> - ujq - ei+k+q + ei) 

XJ\ 2 ^ ( u> B(l,l + q) x 



(19) 



E 

tq 



x 0(- e/ )0(e ;+fc )0(-Q +9 ) 

X 8(u> - UJQ - £l +k + f-l+q), (20) 

where e k = — Jcos(fc) and B(k,l) = sin(fc) — sin(i). 9(x) 
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denotes the step function, and we assume ujq > 2J. In 
this approximation, phonon-assisted particle-hole exci- 
tation starts from the second order of A and is given 
by xi 2a ^(fc,w) and x^c b \k,uj). Thus this excitation is 
expected to increase as the power law X b with 6 = 2. 
This behavior explains the power-law behavior of the in- 
tegrated intensity of the excitation with (AJ/wo) 1 ' 8 as 
shown in the inset of Fig. [Ha). 

Xc 2 °^ and Xc 2& ^ include the particle-hole excitation ac- 
companied by a phonon with the energy of ujq. The 
particle-hole excitation can scan the full-energy range of 
charge excitation with the width of 2 J. Therefore, the 
phonon-induced spin excitation in the XY spin-Peierls 
model is expected to be located between wo and ujo + 2 J, 
which is seen in Fig. U] and Fig. [DJb). At k = n,by us- 
ing eqs. (fTT))) and (|2T)|) . we obtain a cusp- like structure 
as shown in the iset of Fig. [5jb), which is consistent 
with the results shown in in Fig. [5J We note that this 
cusp-like structure is also similar to the phonon-assisted 
magnon absorption observed in the ID Mott insulator- 
s' CUO3. [42|, [43[ Thus we conclude that the full-energy- 
range scanning with particle-hole excitation accompanied 
by a phonon is important to understand the phonon- 
induced spin excitation. 

2. Phonon excitation spectrum 

We examine how phonon excitations arc influenced by 
spin-phonon coupling. Figure Eta) shows the phonon ex- 
citation spectrum for the Einstein phonon with ujq = 3 J. 
The introduction of spin-phonon coupling with XJ/ujq = 
0.5 into a frustrated spin model (ao = 0.36) changes the 
phonon excitation spectrum from Fig. [6ja) to Fig. [6jb). 
Two changes are seen in the phonon excitation spectra. 
One is a slight shift of the phonon main peak toward 
higher energy near q = tt, accompanied by small hump 
at the high energy side as shown in the inset. The en- 
ergy position of the hump structure is the same as that of 
phonon-induced spin excitations seen in Fig. [DJb). The 
other change is the emergence of low-energy phonon com- 
ponents (a; < 1.5J). The strongest change in intensity 
appears near uj = at q ~ tt. We note that, since the 
phonon excitation occurs without spin flipping, i.e., in 
the Hilbert space of zero total spin, the energy positions 
of the spin-induced phonon excitations are different from 
those of spin excitations shown in £HV Bl Both the high- 
and low-energy changes predominantly occur near q = it, 
which is due to the nature of spin-phonon coupling as 
shown in Appendix. We find that these changes are in- 
sensitive to the presence of ao and are general features 
of the spin-phonon coupling. 

V. SUMMARY 

We have investigated spin excitation for a spin-Peierls 
chain with nearest-neighbor and next-ncarest-ncighbor 




2.0 m/J 4.0 6.0 

FIG. 6. (Color online) Phonon excitation spectrum in a 
16-site spin-Peierls chain (ao = 0.36, ujq /J — 3.0) without 
spin-phonon coupling \J/ujo = (a) and with the coupling 
AJ/wo = 0.5 (b). The inset shows the tail of main peak (red 
solid line) together with a single Lorentzian curve obtained 
by fitting the main peak (blue dotted line). The difference 
between the two lines in the inset demonstrates the presence 
of hump structure at the high-energy side of the main peak. 



Heisenberg spin exchange interactions, taking CuGc03 
into consideration. We consider a gapped and dispersion- 
less (Einstein) phonon as the lattice degree of freedom. 
We then apply a dynamical density matrix renormaliza- 
tion group method to calculate dynamical spin correla- 
tion function at zero temperature. 

We have found a new spin excitation assisted by non- 
softening phonon at the energy region higher than the 
phonon energy. The new spin excitation shows a disper- 
sive feature with strong intensity near k = it. There is 
no corresponding structure in the effective J\-Ji model 
derived from the spin-Peierls model as expected. This 
demonstrates the importance of treating the phonons 
quantum- mechanically. 

We have also found the new spin excitation appears in 
the XY spin-Peierls model. This model is equivalent to 
the Su-Schricffcr-Hccgcr model. We have shown that the 
behaviors of the new excitation are explained by charge- 
charge correlation function in the SSH model. 

The phonon excitation spectrum is also influenced by 
the spin-phonon interaction. We have found the shift of 
main phonon peak toward higher energy side near q = tt, 
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accompanied by a new tail structure whose energy range 
is the same as the new spin excitation. In addition, new 
spin-assisted phonon structures appear at low-energy re- 
gion. The fact that the strong modification occurs near 
q = tt is understood by taking into account the momen- 
tum dependence of the spin-phonon coupling. 

The effect of spin-phonon coupling on inelastic neutron 
scattering has been discussed so far in literature. [39| 
However, there is no work demonstrating the intensity 
distributions for the spin-Pcicrls model as far as we know. 
Therefore, we believe that the present results will be help- 
ful for analyzing inelastic neutron scattering data in spin- 
phonon coupled systems. 

In CuGeOs, inelastic neutron scattering experiments 
have clearly revealed the presence of both lowest-energy 
branch of spin excitation with strong intensity and multi- 
spinon continuum. [27i - [29| This is consistent with the re- 
sults of the J1-J2 model. Experimental data has also 
shown high-energy structures above the continuum. If 
the structures are originated from spin degrees of free- 
dom, they may be due to the coupling with phonon. It 
is desired to resolve the structure into phonon and spin 
components by either detailed analyses of momentum de- 
pendence of their intensity or new polarized inelastic neu- 
tron scattering experiments. 
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Appendix A: Momentum dependence of 
spin-phonon coupling 

To examine the nature of the spin-phonon coupling 
term ([J), we introduce Holstcin-Primakoff bosons assum- 
ing two sublatticcs. The corresponding boson operators 
are represented by A m and B m with integer m. We con- 
sider the Fourier transformation of the boson operators: 

S fe = \fh ' ' Jl ' "j witn «2m-i = A m , a 2m = B m . By 

neglecting fourth- or more higher-order terms, the spin- 
phonon term (j4]) is rewritten by 



(Al) 



with 



KM 



2VZ- 



k 

sin(g)a[a fc _ g + sin(/c)(aJ.al fc+9 - a fc a_ fc _ 9 ; 



(,A2) 



where 6^ and b q are the creation and annihilation opera- 
tors of phonons in the momentum representation. 

We can find that the q = phonon has no coupling 
with spin in eq. (|A1|I . Actually, the term including the 
q = 0-phonon operators reads 

K' SP (<1 = 0) = *^=<>o + &o) sin ( k )(a\al k - a fe 5_ fc ) 



AJ 



-n'(g = 0) = 0. 



afe«-fc) 
(A3) 



However, the interaction of the q = it phonon and spin 
remains finite: 

H' sp (q = tt) = + bl) E cos ( fc ) 

X ( S fc+7r/2 S -fc+7r/2 ~ a fc+7T/2a-fc+7r/2) 



*0. 



(A4) 



This is the reason why phonons near q = tt are affected 
strongly by the spin-phonon coupling as compared with 
the region near q = 0. 
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